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SECTION  I 
INTRODUCTION 

SCX  is  a two  dimensional  ground  burst  EMP  (electromagnetic 
pulse)  computer  code.  The  general  numerical  methods  used  in  the 
code  are  documented^  elsewhere  and  will  not  be  discussed  here. 
Briefly,  the  code  obtains  the  solution  to  Maxwell’s  equations  in 
the  source  region  of  i surface  nuclear  burst.  Because  the  solu- 
tion is  obtained  in  the  source  region,  several  nonlinearities 
are  inherent  to  the  problem.  First,  the  conductivity  of  the 
medium  depends  strongly  on  the  total  electric  field.  The  effect 
has  always  been  modeled  in  the  code.  Second,  the  source  terms 
are  themselves  influenced  by  the  fields.  This  effect  is  generally 
referred  to  as  "self-consistency",  and  until  recently  was  not  in- 
cluded in  the  SCX  code.  This  paper  reviews  the  methods  used  to 
•model  self-consistency  in  SCX,  and  presents  comparative  results 
of  calculations  before  and  after  the  effect  was  included  in  SCX. 

An  exact  representation  of  self-consistency  requires  the 
solution  to  the  equations  of  motion  for  the  Compton  electrons . 
Clearly,  for  an  EMP  computer  code  with  two  space  dimensions 
plus  time,  this  is  impractical.  The  amount  of  storage  required 
is  not  available,  and  the  running  time  would  render  the  code 
economically  useless.  Fortunately,  methods  have  been  devised 
which  allow  for  the  inclusion  of  the  self-consistent  effect  in 
an  approximate  fashion s . These  methods  require  a minimal  amount 
of  storage  and  cause  only  slight  increases  in  running  time. 

The  sources  of  the  EMP  are  the  Compton  recoil  electrons 
created  through  the  device  radiation  interactions  with  the  at- 
mosphere. In  SCX  these  sources  are  described  as  current  den- 
sities in  the  radial  and  transverse  directions.  The  current 
densities  were  obtained  from  Monte  Carlo  transport  calculation 
results  which  were  then  curve  fit  for  use  in  the  code.  The 
transport  results,  being  completely  independent  of  the  EMP  cal- 
culation, do  not  contain  any  effects  due  to  interactions  with 
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electromagnetic  fields.  To  include  the  self-consistent  effect, 
some  modification  must  be  made  to  the  source  terms  within  the 
SCX  code.  This  leads  to  several  necessary  approximations,  the 
impact  of  which  will  be  discussed  below. 
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SECTION  II 

A DESCRIPTION  OF  THE  SELF-CONSISTENCY  MODEL 


The  self-consistency  model  used  is  derived  from  EMP  Theo- 
retical Note  77,  Volume  2-4,  by  H.  J.  Longley. ^ The  note 
describes  a way  of  modifying  a purely  radial,  analytic  current 
source  to  obtain  self-consistent  radial  and  transverse  currents. 
The  method  is  based  on  electron  turning  in  the  presence  of 
electromagnetic  fields. 


To  determine  the  amount  of  turning,  a group  of  electrons  is 
followed  in  various  time  constant  electromagnetic  environments. 
These  electrons  are  recoils  created  by  Compton  scattered,  mono- 
energetic  gamma  rays  and  are  chosen  to  represent  a physically  rea- 
listic distribution  of  Compton  recoil  angles  and  energies.  The 
equations  of  motion  for  the  electrons  are  differenced  and  solved 
numerically.  The  computation  proceeds  in  time  until  an  electron’s 


kinetic  energy  is  within  17o  of  its  rest  maSs  energy.  In  STP  air, 

o 

1 MeV  electron  has  a range  of  .49  g/cm  . An  electron  with  a ki- 


netic energy  equal  to  107<>  of  J rest  mass  energy,  has  a range 
2 

of  .0049  g/cm  . So  a 1 MeV  eL  on  slowed  to  17.  of  its  rest 


mass  energy  is  easily  within  .17,  of  its  final  range.  At  this 
point,  the  electron's  final  radial  and  transverse  positions  are 
recorded.  An  average  is  taken  of  the  final  positions  for  the 
group  of  electrons  and  these  averages  are  used  to  obtain  self- 
consistent  currents.  The  validity  of  this  method  depends  on  the 
lifetimes  of  the  electrons  and  the  time  steps  used  in  the  SCX 
calculation.  This  matter  will  be  discussed  later. 


To  obtain  the  self-consistent  radial  current,  the  original 
radial,  analytic  current  is  multiplied  by  DX/Rm^  where  DX  is 
the  average  electron  final  radial  position  and  R^^  is  the  mean 
forward  range  of  the  electron  in  the  absence  of  fields.  The 
self-consistent  transverse  current  is  obtained  by  multiplying 
the  original  radial  current  by  DY/Rm^,  where  DY  is  the  average 
electron  final  transverse  position.  In  the  earth's  magnetic 
field,  the  Larmor  radius  of  a 1 MeV  electron  is  about  100  times 
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its  range.  Since  typical  EMP  fields  produce  much  greater  effects, 
such  as  reversing  the  transverse  current  obtained  from  Monte  Carlo 
transport  calculations,  the  geomagnetic  field  will  be  neglected. 

Two  items  are  important  in  this  method.  The  first  is  the 
initial  kinetic  energy,  Eg,  of  an  electron  to  be  tracked.  This 
energy  depends  on  the  initial  gamma  energy,  E^o,  and  the  scatter- 
ing angle  and  is  obtained  directly  from  the  Klein-Nishina  equation. 
The  second  item  of  importance  is  the  calculation  of  Rmf  Rm£  is 
the  mean  forward  range  obtained  by 


= i-  f 


where 
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and  oe  is  the  angular  differential  cross-section  obtained  from  the 
Klein-Nishina  formula,  dag  is  the  solid  angle  associated  with  the 
scattering  angle  0£  of  the  recoil  electron,  0e  is  the  angle  between 
the  initial  direction  of  propagation  of  the  gamma  and  the  direction 
of  the  electron's  recoil.  R is  the  range  obtained  from  a fit  to  ex- 
perimental mean  range  versus  energy  data.  The  energy  used  to  obtain 

R is  E„  which  is  a function  of  E and  0 . 
e yo  e 

Our  method  of  obtaining  self-consistent  currents  is  different 
in  several  respects  from  the  method  described  in  EMP  Theoretical 
Note  77.  Where  Longley's  method  used  only  an  analytic,  radial 
current  source;  the  current  sources  used  in  SCX  have  both  radial  and 
transverse  components.  The  general  method  described  in  Note  77  is 
designed  for  use  with  an  analytic  current  source.  The  source  terms 
in  SCX  are,  however,  not  analytic,  having  been  obtained  through  curve 
fits  to  the  results  of  gamma  and  neutron  Monte  Carlo  transport  calcu- 
lations. The  source  terms  serve  as  inputs  to  SCX,  and  are  expressed 
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as  total  currents  in  the  radial  and  transverse  directions.  In 
order  to  adapt  the  general  method  to  our  purposes,  tables  similar 
to  Longley’s  were  generated.  However,  our  tables  are  for  electrons 
recoiling  in  the  same  direction  as  the  initial  gamma  propagation 
direction  and  are  not  averages  of  electrons  recoiling  at  different 
0 . This  was  done  because  the  transport  calculations  which  pro- 
vide the  current  sources  for  SCX  already  include  angular  scatter 
effects  and  the  electron  energy  spectrum  is  already  folded  in. 

In  applying  DX/R  and  DY/R  factors  to  the  SCX  transport  derived 
currents,  first  a total  initial  current  is  calculated  from  the  ini- 
tial transport  derived  radial  and  transverse  currents.  This  total 
current  is  then  treated  in  the  same  manner  that  the  analytic  ra- 
dial current  source  is  treated  in  Longley's  method.  To  accomplish 
this,  the  angle  between  the  positive  radial  axis  and  the  total  ini- 
tial current  is  used  to  transform  the  radial  and  transverse  electric 
fields  to  a new  primed  coordinate  system  where  the  total  initial 
current  is  parallel  to  the  primed  positive  radial  axis  (i.e.,  a 
transformation  to  a coordinate  system  in  which  the  transverse  current 
is  zero).  The  DX/R  and  DY/R  factors  are  applied  to  the  total  ini- 
tial current  and  the  resulting  primed  self-consistent  radial  and 
transverse  currents  are  transformed  back  to  the  original  coordinate 
system  to  obtain  the  final  self-consistent  currents. 

In  addition  to  following  single  electrons  rather  than  proba- 
bilistically representative  groups,  our  method  differs  from  Longley's 
in  two  respects.  First,  rather  than  using  the  initial  electron 
kinetic  energy  calculated  directly  from  the  Klein-Nishina  equation, 
a mean  initial  electron  energy  is  used.  This  energy  is  calculated  by 

tt/2 

Ee  - h f Ee°e  dBe 

V° 

Where  Eg  is  a function  of  0g  and  initial  gamma  energy,  E^o,  and  is 
calculated  from  the  Klein-Nishina  equation. 
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Secondly,  the  range  we  use  is  different  from  R the  mean 
forward  range  used  by  Longley.  We  use  a range  R , which  is 
the  electron  range  calculated  by  the  electron  tracking  sub- 
routine with  all  fields  set  to  zero.  This  range  differs  from 
the  range  obtained  by  using  E and  the  R in  Eqn.  (1)  only  be- 
cause of  roundoff  error. 

In  Longley 's  scheme,  tables  of  DX  and  DY  were  generated  for 
various  field  values  and  gamma  energies.  These  tables  were 
then  fit  by  analytic  functions  and  the  functions  used  to  intro- 
duce self-consistency  into  the  LEMP  code.  In  our  case,  tables 
of  DX/Rc  and  DY/R^  are  generated  and  used  directly  by  SCX, 
along  with  some  interpolation  coding,  to  obtain  self-consistency. 
The  interpolation  scheme  is  simple-minded  and  chosen  to  supply 
smooth  sources  to  SCX. 

The  interpolation  is  basically  as  follows.  The  three  field 
values  calculated  by  SCX  are  the  radial  electric  field  (Er> , 
the  transverse  electric  field  (Efi)  and  the  phi  magnetic  field 
(B^).  These  fields  can  be  thought  of  as  the  three  coordinates 
of  a field  space.  For  the  DX/Rc  tables  each  entry  in  the  table 
represents  a point  in  the  field  space.  Similarly  for  DY/R^ . 

For  a given  set  of  field  values  calculated  by  SCX,  the  inter- 
polation coding  determines  what  eight  points,  the  vertices  of 
a rectangular  solid,  corresponding  to  given  DX/Rc  or  DY/RC 
surround  the  point  P whose  coordinates  are  given  by  the  three 
field  values.  This  rectangular  solid  can  be  broken  up  into 
eight  sub-solids  by  passing  three  planes  through  the  fields 
value  point.  The  planes  are  parallel  to  the  six  faces  of  the 
original  solid  and  generate  one  sub-solid  for  each  vertex  of 
the  original  solid.  The  interpolation  scheme  weights  the  DX/Rc 
or  DY/RC  value  at  a particular  vertex  by  a volume  obtained  by 
subtracting  the  volume  of  the  sub-solid  of  that  vertex  from 
the  volume  of  the  original  solid.  This  weighting  is  done  for 
all  eight  points,  summed  and  divided  by  the  total  volume  of 
the  original  solid.  The  scheme  is  smooth  and  has  the  advantage 
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that  if  the  fields  values  fall  exactly  on  a point  in  the  tables, 
the  precise  values  of  DX/Rc  and  DY/Rc  from  the  tables  are  ob- 
tained. 

The  equation  of  motion  used  in  the  electron  tracker  sub- 
routine is: 

p = - 1 e | j£  + v x §j  - A -2— 

1 1 ipi 

where  p is  the  electron's  momentum,  e the  electron's  charge, 
v the  electron's  velocity,  £ the  electrical  intensity  of  the 

■4- 

environment,  B the  magnetic  intensity,  and  A is  a slowing  term 
like  ^ which  includes  energy  losses  due  to  ionization,  multi- 
ple scattering  and  radiation.  If  £ and/or  B are  large  enough 
their  contributions  will  overcome  the  energy  loss  term  A and 
the  electrons  will  never  come  to  rest.  These  are  termed  run- 
away electrons.  The  self-consistency  model  used  here  includes 
a range  of  field  values  which  generate  DX/Rc  and  DY/R£  tables 
that  exclude  run-away  electrons.  Therefore,  the  interpolation 
coding  holds  field  values  to  the  limits  used  in  generating  the 
tables . 

In  estimating  the  effect  of  self-consistency  on  the  con- 
ductivity we  have  to  consider  the  effects  of  the  electric  fields 
on  an  electron's  kinetic  energy,  since  the  kinetic  energy  de- 
termines the  amount  of  ionization.  If  an  electron  has  an  ini- 
tial velocity  in  a given  direction  an  electric  field  parallel 
to  the  velocity  vector  will  increase  or  decrease  the  electron's 
kinetic  energy  depending  on  the  sign  of  the  field. 

Generally,  the  direction  of  the  radial  electric  field  is 
positive,  away  from  the  burst  source  point.  Similarly  for 
the  Compton  recoil  electrons.  Therefore,  the  radial  electric 
field  tends  to  reduce  the  recoil  electron  energy  and  thereby 
reduce  the  ionization  due  to  electrons.  Initially  the  theta 
electron  velocity  in  the  transformed  system  is  zero  so  that  the 
theta  electric  field  will  increase  the  electron's  theta  momentum 
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regardless  of  the  field's  sign.  These  approximate  arguments  lead 
to  the  following  correction  factor,  f , to  the  ionization  rate. 


f = 

q 


E + W 
e 


where 

W = j e | j|^DY|  - ErDxj 

and  |e|  is  the  absolute  value  of  the  electron's  charge.  W is  an 
estimate  of  the  work  done  by  the  fields  on  the  electron  and  there- 
fore changes  the  energy  available  for  ionization. 

The  applicability  of  this  self-consistent  scheme  is  question- 
able when  the  time  steps  used  by  fields  code  differencing  are  com- 
parable or  less  than  the  lifetimes  of  the  electrons.  In  real  time, 

- 8 

for  gammas  of  1.5  MeV,  electron  lifetimes  are  on  the  order  of  10 
second.  In  retarded  time,  due  to  turning,  this  time  may  be  much 
larger  since  there  is  a component  of  the  electron  velocity  which  is 
parallel  to  the  gamma  wave  front. 

In  a typical  SCX  run,  the  time  steps  during  the  prompt  gamma 

-9  -8 

peak  are  10  second.  After  the  peak,  time  steps  are  10  second 

and  larger.  By  examining  electron  trajectories  for  typical  SCX 

environments  it  is  apparent  that  the  electrons  frequently  turn  back 

and  complete  loops.  But  it  is  still  likely  that  the  final  position 

of  an  electron  is  in  the  same  general  direction  from  the  electron's 

original  position  as  the  position  of  an  electron  at  the  end  of  a 

time  step  shorter  than  the  electron's  lifetime.  Since  the  electron 

is  slowing  down,  we  expect  that  the  electric  fields,  at  least,  will 

have  more  effect  on  the  electron's  position  near  the  end  of  its  life 

than  at  the  beginning  where  it  has  large  kinetic  energy.  In  the 

present  SCX  calculations,  the  self-consistent  effect  is  probably 

exaggerated  during  the  prompt  gamma  peak. 


.it 


An  approximate  correction  in  such  situations  might  be  to 
scale  the  turning  by  a factor  tc/t 
and  t 


'el 


0,  where  t is  the  time  step 

S cX.  s 

is  the  electron  lifetime.  A better  factor  would  be 


ft  ^s^eJjP 


The  ffc  fac- 


which  more  heavily  weignts  time  steps  close  to  tg£. 
tor  is  plausible  because  the  non-relativistic  equation  for  a 
displacement  s due  to  a constant  force  on  a mass  m is 


s = 


The  ft  factor  would  scale  the  angle  that  the  position  vector  of 
the  electron's  final  position  makes  with  the  initial  gamma  propa- 
gation direction.  To  accomplish  this,  take  the  original  DX/R 


and  DY/R  . Compute 


a = [DY/R  ] / [DX/R  ] = tan  6 


m 


= i [DX/R  ]2  + [DY/R  ]2  I % 


8 = arctan  ja  • f t j 


[DX/R  ]'  = m cos  (8) 


[DY/Rc] ' = m sin  (8) 

and  use  [DX/R£] ' and  [DY/Rc]'  as  before. 


10 


I 


I 


f 


i 

;> 


4 

i 


3 

r!K 


.£ 


I 


1 1 


% 

r- 


1 


I S* 

t 1 


| 

I 


$ 

% 


SECTION  III 

THE  EFFECT  OF  SELF-CONSISTENCY  ON  THE  CURRENTS 

To  interpret  the  plots  overlaying  self-consistent  and  non- 
self-consistent time  histories  of  SCX  runs,  it  is  useful  to  first 
describe  the  individual  effects  of  the  fields  on  a single  electron. 

-<-  -4r.  -ir 

For  an  electric  field  E,  the  force  on  a charge  q is  F = E q. 

A 

An  electron  in  a positive  Er  field  experiences  a force  in  the  -r 
direction  which  contributes  to  a positive  radial  conventional 
current.  Since  the  Compton  recoil  electrons  are  streaming  radially 
outward,  constituting  a negative  conventional  current,  the  iso- 
lated effect  of  a positive  E field  is  to  reduce  the  magnitude 
of  the  negative  radial  current  -J  . The  same  type  of  argument 
indicates  the  effect  of  EQ  on  the  theta  current. 

SCX  calculates  a which  is  negative.  Since  | J^_  | is  usually 
greater  than  |JQ|,  we  first  consider  the  effect  of  a magnetic 
field  on  a purely  radial  current. 

•f-  ^ 

An  electron  with  a velocity^v  in  the  +r  direction  will  experi- 
ence a magnetic  force  F = q(v  x b).  For  negative  F will  be  in 
the  -0  direction  and  so  contribute  to  a positive  theta  conventional 
current . 

Occasionally,  J^.  and  JQ  are  the  same  order  of  magnitude.  In 
the  extreme  case  where  the  conventional  current  is  purely  in  the  -0 
direction,  the  magnetic  field  will  contribute  to  a negative  radial 
conventional  current. 

In  determining  whether  it  is  an  electric  field  or  the  magnetic 
field  which  dominates  the  electron  turning,  it  is  useful  to  be  able 
to  make  rough  comparisons  between  the  effects  of  the  electric  and 
magnetic  fields. 

For  an  electron  of  initial  energy  Ee>  initial  speed  v , and 
absolute  charge  e,  we  define  Fm  to  be  the  maximum  magnetic  force 
on  the  electron  and  F to  be  the  maximum  slowing  force  on  the 
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Since  these  two  forces  are  mono toni cal ly  increasing 
functions  of  the  electron  velocity,  (except  for  F when  the  elec- 

b 

tron  energy  is  below  0.5  MeV)  the  maxima  occur  at  the  electron’s 


maximum  speed,  i.e. 

.75  MeV  and  v is  2.75  x 10 
o 

pared  as  follows: 


For  an  incident  gamma  of  1.5  MeV,  E is 


8 


m/sec . 


F„  and  F can  then  be  com- 
m s 


F /e  = B.  • 2.75  x 10c 

m <p 


A(v) 


- = 3.49  x 105 


F /e  = 

s e 


The  slowing  term  A(vq)  is  obtained  from  a fit  to  experimental 


data  of  electron  energy  as  a function  of  electron  mean  range. 

dE 

fit  is  differentiated  with  respect  to  range  to  obtain  . 


This 


With  the  use  of  F and  F„  we  can  predict  the  combination  of 

ms  r 

effects  of  the  various  fields  on  the  currents. 

The  radial  current  overlays  in  Figs.  3 and  12  show  that  for 
ranges  of  500m  and  1000m  the  self-consistent  model  reduces  the  mag- 
nitude of  the  radial  current  until  past  a microsecond.  This  is 
especially  noticeable  at  the  time  of  the  prompt  gamma  peak  and  also 
past  10  shakes  where  the  reduction  in  J increases  markedly  with 
time  until  near  a microsecond.  The  2000m  radial  currents  overlay 
exactly,  which  prompts  us  to  consider  the  close-in  ranges  and  the 
2000m  range  separately.  Evidently,  at  2000m  for  this  yield,  the 
fields  are  reduced  enough  to  show  only  small  self-consistent  effects, 
primarily  in  J0  and  E0. 

As  mentioned  earlier,  the  effect  of  a positive  Er  field  is  to 
reduce  the  magnitude  of  a negative  Jr.  For  the  two  close-in  ranges, 
Er  is  positive  throughout  the  calculation.  Also,  a negative  B^  will 
only  increase  the  magnitude  of  negative  Jr  when  J0  is  negative  and 


J 


9 


In  the  self-consistent  case,  we  see  that  |Jr|/Uel  is 


close  to  unity  at  usee.  At  this  point,  in  Fig.  3,  the  self- 
consistent  Jr  is  increased  as  expected.  By  comparing  the 
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I 
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'm 


-I 
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non-self-consistent  JQ  and  Jr,  it  is  seen  that  separation  between 
the  self-consistent  J and  non-self-consistent  Jr  continues  to 
increase . 

The  more  pronounced  -J  reduction  at  the  time  of  the  prompt  gamma 
peak  is  due  to  the  combined  peaking  of  E^  and  B^.  After  the  prompt 
gamma  pulse,  whereas  Er  is  saturated  and  remains  reasonably  constant 

out  to  neutron  arrival  at  around  10  microseconds,  B.  steadily  increases 

<P 

in  magnitude  and  thus  increases  the  separation  of  the  Jr  overlays. 

It  is  interesting  to  note  that  at  500m  the  self-consistent  B^  starts 
leveling  off  at  around  2 usee  and  then  starts  to  decrease  at  about 
5 psec.  The  separation  in  the  Jr  overlays  follows  this  behavior 
until  the  non-self-consistent  J0  becomes  larger  than  the  non-self- 
consistent  J . Similar  behavior  is  shown  at  1000m. 

At  2000m  there  is  no  visible  effect  of  the  self-consistent 

model  on  the  radial  current.  This  is  plausible  on  the  basis  of 

rough  field  comparisons.  The  maximum  value  of  E at  this  range 

is  5 x 10  . The  maximum  absolute  value  of  B.  is  3 x 10_£+. 

4 9 5 

In  this  case,  Fm/e  is  - 8 x 10  . Since  F g/e  is  ~ 3.5  x 10  , it  seems 

reasonable  that  the  self-consistent  effect  on  J due  to  E will 

r r 

be  negligible,  and  the  effect  due  to  B^  will  be  small,  particularly 
since  at  this  range,  | JT^.  I > > I J0 1 and  nearly  all  the  kinetic  energy 
of  the  electron  is  in  the  radial  direction. 

The  self-consistent  theta  currents  for  the  two  close-in  ranges 
show  three  interesting  features.  First,  while  the  non-self-consistent 
theta  currents  are  always  negative,  the  self-consistent  theta  currents 
are  nearly  always  positive.  Second,  the  self-consistent  theta  currents 
follow  the  prompt  gamma  pulse  in  a much  more  obvious  fashion  than  the 
non-self-consistent  theta  currents.  Third,  after  the  prompt  gamma 
pulse,  the  self-consistent  theta  currents  dip  and  then  exhibit  a 
gentle  bump,  and  finally  change  sign  after  10  microseconds. 

The  self-consistent  theta  currents  are  nearly  always  positive 
because  B^  nearly  always  predominates  over  E^  and  the  non-self- 
consistent  J is  nearly  always  greater  than  the  non-self-consistent 


m 


m 


Jg.  There  is  a very  short  span  of  time  at  very  early  times  where 
Eg  predominates.  If  values  of  E0  and  F /e  are  compared  at  3 shakes 
in  the  usual  manner,  the  E0  dominance  can  be  shown.  In  this  tiny 
region  of  E0  dominance,  the  self-consistent  theta  currents  are 
negative.  This  situation  is  shown  in  Figs.  5 and  14. 

That  the  theta  currents  are  almost  entirely  determined  by 
is  further  demonstrated  by  the  jagged  time  behavior  of  E0  and  con- 
sequent smearing  of  the  prompt  gamma  pulse.  In  contrast,  the  theta 
currents  are  smooth  and  follow  the  gamma  pulse  quite  well  because 
close-in  the  shape  of  the  gamma  pulse  is  preserved  in  B^. 

Finally,  the  dip  and  gentle  bump  behavior  is  exhibited  in  B^ 
but,  due  to  the  obvious  non-linear  relationship  of  electron  turning 
to  the  magnitude  of  B^  , the  similarity  of  shape  between  the  self- 
consistent  J0  and  B^  is  not  compelling,  especially  as  the  waveforms 
approach  10  ysec  where  the  non-self-consistent  JQ  becomes  comparable 
to  or  greater  than  the  non-self-consistent  J . Beyond  10  ysec  J0 
crosses  over  due  to  the  fact  that  the  non-self-consistent  J_  becomes 
comparable  to  or  greater  than  (at  500m)  the  non-self-consistent  J . 

In  this  region  the  effect  of  B^  is  to  increase  - Jr,  as  explained 
above,  and  so  EQ  dominates  JQ  behavior.  E0  starts  its  dominance 
before  10  ysec.  The  effect  is  to  reverse  J0 . At  12  ysec  and  500m 
(Fig.  9)  and  22  ysec  and  1000m  (Fig.  18)  E0  crosses  over  and  becomes 
positive.  This  causes  the  self-consistent  J0  to  hump  over  as  it 
heads  for  another  cross-over. 

At  2000m,  rough  field  comparison  shows  that  E0  should  dictate 
J0  behavior.  First,  a vestige  of  the  gamma  pulse  is  seen  in  J0 . 

The  shape  of  the  gamma  pulse  is  preserved  in  E0  but  not  in  B^. 

Second,  J0  is  uniformly  negative  as  is  E0  past  the  start  of  the  gamma 
pulse.  Later  in  time  B^  rises  faster  than  E0  and  at  its  peak  there  is 
a corresponding  dip  in  J0  because,  as  shown  earlier,  a negative  B^ 
acting  on  a negative  Jr  contributes  to  a positive  JQ . 

Now  it  remains  to  examine  the  effects  of  the  self-consistent 
model  on  the  fields  and  the  conductivity. 
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SECTION  IV 

THE  EFFECT  OF  SELF-CONSISTENCY  ON  THE  FIELDS 

A few  words  should  be  said  regarding  the  occasional  ragged- 
ness of  some  of  the  fields.  By  examining  range  plots  of  Eg  , 
it  is  clear  that  choosing  the  inner  boundary  condition  EQ  = 0 
is  inappropriate.  Originally  this  condition  was  chosen  with 
the  assumption  that  the  inner  boundary  is  a perfect  conductor. 
This  assumption  is  certainly  inconsistent  with  the  use  of  non- 
zero theta  currents  at  the  inner  boundary.  Range  plots  of  EQ 
at  early  times  show  a drastic  discontinuity  between  the  inner 
boundary  and  the  first  point  out  in  range.  In  fact,  EQ  is 
increasing  in  an  exponential  fashion  toward  the  inner  boundary 
rather  than  decreasing  to  zero.  After  a few  time  steps,  this 
discontinuity  develops  into  oscillations  of  E0  in  range.  In 
turn  these  oscillations  affect  JQ  which  feeds  back  into  Eg. 

To  minimize  these  oscillations,  a range  current  smoother  has 
been  installed  in  SCX.  This  stopgap  measure  is  helpful  but 
not  completely  effective  as  can  be  seen  in  the  time  plots  o' 
near  the  prompt  gamma  pulse. 

In  addition,  the  calculation  of  the  conductivity  involves 
using  a field  dependent  electron  mobility  which  is  clearly 
affected  by  the  erratic  behavior  of  Eg.  The  conductivity's 
slightly  ragged  behavior  is  fed  back  into  Er  and  into  B^. 

This  problem  should  be  cleared  up,  if  not  eliminated,  by  a 
more  physically  realistic  choice  of  inner  boundary  condition 
for  Eg,  possibly  something  as  simple-minded  as: 


Ee  = "Je/a 


The  three  equations  of  importance  in  SCX  are,  in  retarded 


time,  at  0 = 90°  (on  the  ground) 
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The  usual  arguments  used  to  predict  the  time  behavior  of  Er 

from  1)  are  as  follows.  With  the  magnetic  term  negligible,  at 

very  early  times  J^,  c and  Er  are  very  small  and  so  a Er  is 

negligible  relative  to  J^..  Hence,  the  Er  behavior  is  predicted 

to  be:  E = --  f J dx. 
r cj  r 


'I 

•t 
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After  a time,  oEr  becomes  comparable  to  -Jr-  Physically,  this 
is  described  as  occurring  when  the  Compton  current  is  cancelled  by 
the  conduction  current.  When  this  condition  occurs  -zp^-  is  negli- 
gible, assuming  the  effects  of  the  airground  asymmetry  have  not 

yet  allowed  B,  to  diffuse  into  the  region  of  interest.  If  J and 

* cjt*  3 Er 

a rise  initially  as  e , E_.  saturates,  i.e.,  — - is  small,  for 

3t 

a > aeQ,  where  eQ  is  the  free  space  permitivity.  For  an  a of 
2 x 10^,  saturation  occurs  where  o > 1.77  x 10  ^ and  the  time  of 
saturation  can  be  determined  by  examining  Figs.  6,  15  and  23.  For 
500  and  1000  meters,  saturation  occurs  before  the  prompt  gamma 
peak.  At  2000  meters,  saturation  never  occurs. 

How  accurately  follows  -J  fo  is  estimated  by  a "relaxation 
time"  which  amounts  to  eQ/a.  If  o is  large  enough,  the  relaxation 
time  is  so  short  that  E^  does  in  fact  follow  J^/a,  most  visibly  at 
late  times  where  Jr/o  changes.  At  far  ranges  or  closer  in  at  very 
late  times,  a is  so  small  that  the  relaxation  time  is  too  large  to 
follow  Jr/o. 


17 


Furthermore,  at  the  close-in  ranges,  500  and  1000m,  satura- 
tion occurs  before  the  prompt  gamma  peak  so  that  Er  also  peaks, 
and  at  far  ranges,  2000m,  saturation  occurs  after  the  prompt  gamma 

peak  so  that  the  peak  is  not  preserved  in  E by  following  J /a, 

—If  ^ ^ 

but  rather  from  E = — — / J dt.  For  E , this  results  in  a peak 

more  broad  and  delayed  in  time  from  the  prompt  gamma  peak.  In 

certain  time  domains,  some  of  these  arguments  are  equally  applicable 

to  E0 . 

At  500m  and  1000m,  saturation  occurs  before  the  prompt  gamma 
peak.  The  plots  of  a at  these  ranges  show  that  the  self-consistent 
model  doesn't  greatly  change  the  conductivity.  However,  E0 , 
through  the  field-dependent  electron  mobility,  introduces  some  small 
jaggedness  into  a. 

Since  Er  = -J^/a  until  past  neutron  arrival  where  a is  greatly 
reduced,  thus  increasing  the  relaxation  time,  it  is  reasonable  that 
the  1/a  dependence  of  Er  greatly  exaggerates  the  jaggedness  in  a. 
Since  Jr  is  reduced  by  self-consistency,  Er  is  reduced  as  well. 

An  interesting  portion  of  the  E^  curve  is  at  and  past  neutron 
arrival.  The  a curves  show  a sharp  discontinuity  in  slope  at  neu- 
tron arrival  and  a subsequent  characteristic  hump.  E exhibits  this 
same  slope  discontinuity  and  an  inverted  hump  out  to  about  30  psec. 

Close  in,  before  the  prompt  gamma  peak,  E„  is  driven  by  the 
3 ^ 

-^(rB^)  term.  At  500m  and  1000m,  this  can  be  readily  seen.  In  both 
self-consistent  and  non-self-consistent  plots  of  E0  there  is  a very 
smooth,  sharp  negative  pulse  which  peaks  at  about  5 shakes.  If 
this  pulse  were  due  to  -J0/o  the  E0  pulse  caused  by  a self-consistent 
J0  would  be  opposite  in  sign  to  the  E0  pulse  generated  by  a non-self- 

consistent  JQ.  Examination  of  the  slope  of  the  B,  curve  shows  that 

0 9 r <p 

the  E0  pulse  is  in  fact  driven  by  ^-(rB^).  Between  3 and  5 shakes, 


9t 


B^  rises  rapidly  and  smoothly  to  a peak.  Since  BA  is  negative,  an 


increasing  ^(rB^)  should  give  a negative  EA  value.  It  is  clear 

'6 


0 ^ \ X.  / 0 L lV/  U 

from  the  non-self-consistent  plots  of  Efi  that  the  pulse  ends  and  a 
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sign  change  occurs  at  the  point  where  B peaks  and  turns  over. 
After  peaks,  the  slope  of  doesn't  dr  anything  of  great 
interest  until  neutron  arrival.  In  the  intervening  interval,  E0 
is  driven  by  -Jq/cs  as  can  be  accurately  verified  by  comparing 
-Jq/o  with  actual  val  s of  E0. 

An  important  difference  shown  in  the  overlay  plots  of  EQ 
at  500m  and  1000m  is  that,  whereas  EQ  remains  positive  for  a 
long  time  after  the  negative  pulse  for  the  non-self-consistent 
case,  Eg  remains  negative  for  the  self-consistent  case.  Here 
E0  is  just  following  -JQ/a. 

At  2000m,  an  interesting  feature  is  that  while  self-consis- 
tency reduces  the  magnitude  of  JQ  due  to  B^,  the  self-consistent 
Eg  is  actually  larger  than  the  non-self-consistent  EQ  at  times 
greater  than  10  shakes.  Here  EQ  is  not  driven  by  JQ. 
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SECTION  V 

RADIATION  ENHANCED  GROUND  CONDUCTIVITY 

In  che  past,  SCX  calculations  have  always  assumed  a uniform 
homogeneous  ground  with  constant  conductivity.  However,  in  the 
real  physical  case,  the  deposition  of  radiation  within  the  ground 
results  in  ionization  which  alters  the  conductivity  from  its 
ambient  value.  The  time  variation  of  the  source  and  the  nature 
of  the  deposition  make  the  ground  conductivity  a function  of  both 
space  and  time.  It  is  important  to  determine  the  effects  of  these 
possible  variations  on  the  EMP  environments  calculated  with  the 
SCX  code. 

To  enhance  calculational  speed  and  efficiency,  the  coding 
in  the  SCX  fields  calculation  has  always  implicitly  assumed  a 
constant  ground  conductivity.  With  the  assumption  of  a variable 
conductivity,  the  differenced  form  of  the  radial  equation  in  the 
ground  becomes ^ 


1 hk  fit\  P hk  / Z(p,z)  fit  U k 


O • • H 7* I £j 

z xj  Mg  IT)  p,, 
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After  the  standard  definition  of  constants,  the  following  result 
is  obtained 


Al.E  hk  - A21.B.k  - A22.B  k = A3.. 
3 Py  J *tj  J *ij-l  3 


However,  the  following  revised  values  of  several  constants  must 
be  used: 


at  _ 1 , _hk  6t 

A1j  ' ? °ij  “g  Z" 


A,  /I  hk-1  6 1 \ - hk- 1 
3 ,3  8r  py 


Z(p,z)  , 6 t \R  k-1  / Z(p,z)  6 t \R  k-1 

2 ~ \ 2 25i7/*ii-l 


y/ 


The  other  equation  to  be  differenced  which  involves  the  con- 
ductivity is 


3B. 

R(p,z)  ^ 


13  1 

F Ip^V  + OUEz  ly  ^ 

g 


Assuming  a variable  conductivity,  the  following  difference  equation 
results 


mim 


R(p,z)  - 


pi  «t\  v k . 1 1 , k «t\F  k 

\.  +l^+WsT-j\. 


1 k-l  6t\T7  k-1  , /D/  N 
7_aij  pg  r K..  + / R(p,z> 


5t\g  k-l 

Pi  + Pi+1  «P) 


' Pi-1  <St\B  k , / pi+l  6t\B  k 

p.  + p.^  6pl  ( Pi  + Pi_i  ip)  4»i+lj 


Definition  of  constants  results  in 


A4.B  k + A5.E  k = A6.  , 
J *ij  J ^ J 


However,  the  following  constants  require  new  definition 

A5 . = ^7  + crk.  p 2"^  A1i  for  this  case) 

J Q S J 


A6.  = 
J 


1 k-l  5 ti,,  k-l  , N 

' ij  “g  Z-K,  R(P'Z) 


6t\„  k-l 


pi  + pi+i  *7  *ij 


, / pi+l  5t\B  k-l  / pi-l  «t\B  k 

i Pi  + p1+1  «pJ  <Di+lj  "l  Pt  + Pi_i  «P/  ♦i-.ij 


When  these  changes  are  included  in  the  field  calculation  subroutine 
of  SCX,  the  effect  of  radiation  enhanced  ground  conductivity  may 
be  examined. 
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Numerous  models  have  been  proposed  to  approximate  the  behavior 
of  the  ground  conductivity  with  dose.  To  estimate  the  nature  of 
the  effect  in  SCX,  it  is  convenient  to  use  a simple  model  suggested 
by  Graham  and  used  by  Jones^  • In  this  approximation 


Q,z)  = o (constant)  + — 1 - X- r-n  Qe 

s s 8.081  x 10iU  10° 


20z 


where 


Q is  the  ionization  rate  at  the  ground 

8.081  x 10^  = 1 rad  air 

in 

a is  the  depth  in  meters  (a  negative  number) , and 


o (constant)  is  the  normal  ground  conductivity. 

s 


For  a source  on  the  ground,  as  is  the  case  in  SCX,  the  deposition 
beneath  the  surface  is  rather  small  for  ranges  greater  than  a few 
hundred  meters.  A typical  value  for  the  ground  conductivity  in  SCX 
calculations  is  0.01  mhos /meter.  Figure  1 shows  the  radiation  en- 
hanced conductivity  as  a function  of  time  5 cm  below  the  surface 
at  a range  of  2.50  meters  in  a typical  SCX  run.  It  can  be  seen 
that  the  values  change  by  at  most  about  50%  near  the  peak.  Near 
the  prompt  peak,  values  of  the  transverse  electric  field  on  the 
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ATTN:  Library 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Aerojet  Electro-Systems  Co.  Div. 

Aerojet- General  Corporation 
ATTN:  Technic  il  Library 
ATTN:  Thomas  I).  Banscoine 

Aeronutrome  Ford  Corporation 
Aerospace  & Communications  Ops. 

ATTN:  Ken  C.  Attinger 
ATTN:  E.  R.  Poncelet,  Jr. 

ATTN:  L.  B.  Linder 

ATTN:  Technical  Information  Section 

Aeronutronic  Ford  Corporation 
Western  Development  Laboratories  Division 
ATTN-  Samuel  R.  Crawford,  MS  531 
ATTN:  J.  T.  Mattingly,  MS  X-22 
ATTN  Library- 

Aerospace  Corporation 
ATTN:  Bal  Krishan 
ATTN:  Melvin  a.  Bernstein 
ATTN:  S.  P.  Bower 
ATTN  Julian  Reinheimer 
A’lTN:  Irving  M.  Garfunkel 
ATTN:  Dr.  B.  Barry 
ATTN:  Norman  D.  Stockwcll 
ATTN:  Library 


Avco  Research  & Systems  Group 

ATTN:  Research  Library  AS30,  Rm. 
ATTN:  W.  Broiling 

Battellc  Memorial  Institute 
ATTN:  Technical  Library 
ATTN:  David  A.  Dingee 
ATTN:  Dr.  L.  E.  Hulbcrt 
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DEPARTMENT  OF  DEFENSE  CONTRAC  1'OIiS  (Continued) 

The  BDM  Corporation 

ATTN:  Technical  Library 

The  BDM  Cor]>oration 

ATTN:  Technical  Library 
ATTN:  Robert  Ii.  Buchanan 
ATTN:  T.  il.  Neighbors 

Bell  Aerospace  Company 
Division  of  Textron,  Inc. 

ATTN:  Carl  B.  Schoch,  Wpns.  Effects  Grp. 

ATTN:  Martin  A.  Henry 
ATTN:  Technical  Library 

The  Bcndix  Corporation 
Communication  Division 

ATTN:  Document  Control 

The  Bcndix  Corporation 
Research  Laboratories  Division 
ATTN:  Technical  Library 
ATTN:  Donald  J.  Niehaus,  Mgr.  Prgm.  Dev. 

The  Bcndix  Corporation 
Navigation  & Control  Division 
ATTN:  Technical  Library 

The  Boeing  Company 

ATTN:  Datid  Kcinle 

ATTN:  David  L.  Dye,  MS  87-75 

ATTN:  Howard  W.  Wicklein,  MS  17-11 

ATTN:  D.  E.  Isbell 

ATTN:  Dr.  B.  Lampriere 

ATTN:  Aerospace  Library 

Booz- Allen  & Hamilton,  Inc. 

ATTN:  Raymond  J.  Chrisncr 
ATTN:  Technical  Library 

Brown  Engineering  Company,  Inc. 

ATTN.  David  L.  Lambert,  MS  18 
ATTN:  John  M.  McSwain,  MS  18 
ATTN:  Technical  Library,  P.  Shelton,  MS  12 

Burroughs  Corporation 
Federal  & Special  Systems  Group 
ATTN:  Angelo  J.  Mauricllo 
ATTN:  Technical  Library 

Calspan  Corporation 

ATTN:  Technical  Library 

Charles  Stark  Draper  Lalioratory,  Inc. 

ATTN:  Technical  Library 
ATTN:  Kenneth  Fcrtig 

Cincinnati  Electronics  Corporation 
ATTN:  Technical  Library 

Computer  Sciences  Corporation 
ATTN:  Technical  Library 

Computer  Sciences  Corporation 
ATTN:  Alvin  T.  Schiff 

Cutler- Hammer,  Inc. 

AIL  Division 

ATTN:  Anne  Anthony,  Central  Technical  Files 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

University  of  Denver 
Colorado  Seminary 

ATTN:  Security  Officer  for  Ron  W.  Buchanon 
ATTN:  Security  Officer  for  Technical  Library 
ATTN:  Security  Officer  for  Fred  P.  Vcnditti 

The  Dikewood  Corporation 
ATTN:  L.  Wayne  Davis 
ATTN:  Technical  Library 

E-Systems,  Inc. 

Greenville  Division 

ATTN:  Library  8-50100 

Effects  Technology,  Inc. 

ATTN:  Edward  John  Steele 
ATTN:  B.  Wcngler 
ATTN:  Technical  Library 

EG&G,  Inc. 

Albuquerque  Division 

ATTN:  Technical  Library 

ESL,  Inc. 

ATTN:  Technical  Library 
ATTN:  William  Metzer 

Experimental  £•  Mathematical  Physics  Consultants 
ATTN:  Thomas  M.  .Ionian 

Fairchild  Camera  & Instrument  Corporation 

ATTN:  Security  Department  for  Technical  Library 

Fairchild  Industries,  Inc. 

Sherman  Fairchild  Technology  Center 
ATTN:  Leonard  J.  Schreiber 
ATTN:  Technical  Library 

The  Franklin  Institute 

ATTN:  Ramie  II.  Thompson 
ATTN:  Technical  Library 

Garrett  Corporation 

ATTN:  Technical  Library 

General  Dynamics  Corporation 
Pomona  Operation 

ATTN:  Technical  Library 

General  Dynamics  Corporation 
Electronics  Division 

ATTN:  Technical  Library 

General  Electric  Company 
Space  Division 

ATTN:  Technical  Information  Center 

ATTN:  Larry  I.  Chasen 

ATTN:  Joseph  C.  Pedcn,  CCF-8101 

ATTN:  Dante  M.  Tasca 

ATTN:  James  P.  Spratt 

ATTN:  Daniel  Edclman 

ATTN:  J.  Ilannabcck 

General  Electric  Company 
Re-Entry  & Environmental  Systems  Division 
ATTN:  John  W.  Palchcfsky,  Jr. 

ATTN:  Technical  Library 


DEPARTMENT  OF  DKKENSK  CONTRACTORS  (Continued) 

General  Electric  Company 

Ordnance  Systems 

ATTN:  Joseph  J.  Rcidl 

General  Electric  Company 

TEMPO-Center  for  Advanced  Studies 
ATTN:  Hoyden  11.  Rutherford 
ATTN:  DASIAC 

General  Electric  Company 

ATTN:  Richard  C.  Fries,  CSP  G-7 
ATTN:  Technical  Library 

General  Electric  Company 

Aircraft  Engine  Group 

ATTN:  John  A.  Ellerhorsl,  E-2 
ATTN:  Technical  Library 

General  Electric  Company 

Aerospace  Electronics  Systems 

ATTN:  George  Francis,  Drop  233 
ATTN:  Charles  M.  Hewison,  Drop  G2-1 
ATTN : Technical  Library 

General  Electric  Company 
ATTN:  Technical  Library 

General  Research  Corporation 
ATTN:  John  Ise,  Jr. 

ATTN:  Robert  D.  Ilill 

ATTN:  Technical  Information  Office 

Grumman  Aerospace  Corporation 

ATTN:  Jerry  Rogers,  Department  533 
ATTN:  Technical  Library 

GTE  Sylvania,  Inc. 

Electronics  Systems  Group-Eastern  Division 
ATTN:  James  A.  Waldon 
ATTN:  Charles  A.  Thornhill,  Librarian 
ATTN:  Leonard  L.  Hlaisdcll 

GTE  Sylvania,  Inc. 

ATTN:  Herbert  A.  Ullman 

ATTN:  Mario  A.  Nurcfora,  11  4 V Group 

ATTN:  S.  E.  Perlman,  A.S.M.  Department 

ATTN:  David  P.  Flood 

ATTN:  Emil  P.  Motchok,  Comm.  Syst. 

Harris  Corporation 

Harris  Semiconductor  Division 

ATTN:  Wayne  E.  Abare,  MS  1G-111 
ATTN:  Carl  F.  Davis,  MS  17-220 
ATTN:  T.  L.  Clark,  MS  1040 
ATTN:  Technical  Library 

Hazcltinc  Corporation 

ATTN:  M.  Waite,  Technical  Information  Center 

Hercules,  Incorporated  * 

ATTN:  Technical  Library 
ATTN:  R.  Woodruff,  100K-20-W 

Honeywell  Incorporated 

Government  4 Aeronautical  Products  Division 
ATTN:  Ronald  R.  Johnson,  A- 1022 
ATTN:  Technical  Library 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Honeywell  Incorporated 

Aerospace  Division 

ATTN:  Stacey  H.  Graff,  MS  725-J 
ATTN:  Technical  Library 

Honeywell  Incorporated 

Radiation  Center 

ATTN:  Technical  Library 

Hughes  Aircraft  Company 

ATTN:  Billy  W.  Campbell,  MSC-E-110 
ATTN:  Technical  Library 

Hughes  Aircraft  Company 

Ground  Systems  Group 

ATTN:  Library,  MS  C-222 

Hughes  Aircraft  Company 

Space  Systems  Division 

ATTN:  Edward  C.  Smith,  MSA-620 
ATTN:  William  W.  Scott,  MS  A-1080 
ATTN:  Technical  Library 

IBM  Corporation 

ATTN:  Frank  Frankovsky 
ATTN:  Technical  Library 

in'  Research  Institute 
ATTN:  ACOAT 
ATTN:  Technical  Library 

HT  Research  Institute 

ATTN:  Irving  N.  Mindel 
ATTN:  Technical  Library 

Institute  for  Defense  Analyses 

ATTN:  IDA  Librarian,  Ruth  S.  Smith 

Intclcon;  Rad  Tech 

ATTN-  it.  L.  Mertz 
ATTN:  Ralph  H.  Stahl 
ATTN:  Technical  Library 

International  Telephone  4 Telegraph  Corporation 
ATTN:  J.  Gulack,  Def.  Sp.  Gip. 

ATTN:  Tcchnical  Library 

Ion  Physics  Corporation 

ATTN:  Ronert  D . Evan3 
ATTN:  II.  Mildc 
ATTN:  B.  Evans 
ATTN:  Technical  Library 

Johns  Hopkins  University 

Applied  Physics  Laboratory 
ATTN:  Technical  Library 

Litton  Systems,  Inc. 

Data  Systems  Division 

ATTN:  Technical  Library 

Litton  Systems,  Inc. 

Guidance  4 Control  Systems  Division 
ATTN:  John  P.  Rctzlcr 
A.TN:  ValJ.  Ashby,  MS  67 
ATTN:  Technical  Library 

Ktech  Corporation 

ATTN:  Dr.  D.  Keller 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Commut'd I 

Hainan  Sciences  Corporation 
ATTN:  W.  Foster  Rich 
ATTN:  Walter  E.  Ware 
ATTN-  John  It.  Hoffman 
ATTN.  Donald  II.  Bryce 
ATTN:  Albert  P.  Bridges 
ATTN:  Frank  H.  Shelton 
ATTN:  T.  Meagher 
ATTN:  Dr.  I).  C.  Sachs 
ATTN:  J.  C.  Nickell 
ATTN:  J.  Oscarson 
ATTN:  D.  Williams 
ATTN:  It.  McClellan 
ATTN:  E.  Walsh 
ATTN:  Dr.  P.  Wicsolmann 
ATTN:  Library 

Litton  Systems,  Inc. 

AMECOM  Division 

ATTN:  Technical  Library 

Lockheed  Missiles  &•  Space  Co.  inc. 

ATTN:  George  F.  Heath,  Dept.  81- M 
ATTN:  Kevin  McCarthy,  0-S5-71 
ATTN:  Hans  L.  Schneemann,  Dept.  81-GI 
ATTN:  L-365,  Dept.  81-20 
ATTN:  Philip  J.  Hart,  Dept.  81-11 
ATI'S:  Benjamin  T.  Kimura,  Dept.  81-U 
ATTN:  D.  M.  Tcllcp,  Dept.  81-01 
ATTN:  Dr.  M.  Miller 
ATTN:  A.  O.  Burford 
ATTN:  Technical  Library 

Lockheed  Missiles  4-  Sp„cc  Company 
ATTN:  P.  G.  Underwood 
ATTN:  Technical  Information  Center,  D/Coll. 

LTV  Aerospace  Corporation 

Vought  Systems  Division 

ATTN:  Technical  Data  Center 

LTV'  Aerospace  Corporation 

Michigan  Division 

ATTN:  James  F.  Sanson,  B-2 
ATTN:  Technical  Library 

M.LT.  Lincoln  Laboratory 

ATTN:  Leona  Loughlm,  Librarian  A-082 

Martin  Marietta  Aerospace 

Orlando  Division 

ATTN:  MonaC.  Griffith,  Lib.,  MP-30 
ATTN:  JackM.  Ashford,  MP-537 

Martin  Marietta  Corporation 

Denver  Division 

ATTN:  Paul  G.  Kase,  Mail  8203 

ATTN:  Ben  T.  Graham,  MS  PO-15J 

ATTN:  Jay  R.  McKee,  Research  Library  6G17 

Maxwell  Laboratories,  Inc. 

ATTN:  Richard  A.  Fitch 
ATTN:  Victor  Fargo 
ATTN:  Technical  Library 

McDonnell  Douglas  Corporation 
ATTN:  Tom  Endcr 
ATTN:  Technical  Library 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

McDonnell  Douglas  Corporation 
ATTN:  W.  R.  Spark,  MS  13-3 
ATTN:  A.  P.  Venditt,  MS  11-1 
ATTN:  Stanley  Schneider 
ATTN:  Dr.  R.  J.  Reck 
ATTN:  Dr.  11.  M.  Berko  wit  7. 

ATTN:  Technical  Library  Services 

McDonnell  Douglas  Corporation 
ATTN:  Thomas  J.  Ltmdregan 

Mission  Research  Corporation 
ATTN:  William  C.  Hart 
ATTN:  Conrad  L.  Longmire 
ATTN:  Daniel  F.  Higgins 
ATTN:  Technical  Library- 

Mission  Research  Corporation 
ATTN:  David  E.  Merewcther 
ATTN:  Larry  D.  Scott 
ATTN:  Technical  Library 

The  '.itre  Coriioration 

ATTN:  Theodore  Jarvis 
ATTN:  Louis  Brickmore 
ATTN:  M.  E.  Fitzgerald 
ATTN:  Library- 

National  Academy  of  Sciences 

ATTN:  R.  S.  Shane,  Nat.  Materials  Advsy. 

Northrop  Corporation 

ATTN:  George  II.  Towner 
ATTN:  Vincent  R.  DcMartino 
ATTN:  John  M.  Reynolds 
ATTN:  Technical  Library 

Northrop  Corporation 

ATTN:  Technical  Library 

Northrop  Corporation 
ATTN:  Library 
ATTN:  David  M.  Pocock 

Palisades  Inst,  for  Research  Services  Inc. 

ATTN:  Records  Supervisor 

Perkm-Elmcr  Corporation 
ATTN:  Technical  Library 

Philco-Aeronut 

ATTN:  R.  Lyons 

Physics  International  Company 

ATTN:  Document  Control  for  Bernard  H.  Bernstein 
ATTN:  Document  Control  for  John  II.  Huntington 
ATTN:  Document  Control  for  Technical  Library 
ATTN:  Dr.  J.  Shea 
ATTN:  K.  Childers 

Proccdync  Corporation 

ATTN:  Peter  Horowitz 
ATTN:  Technical  Library 

Prototype  Dev.  Asso. 

ATTN:  T.  McKinley 

R 4-  D Assoc. 

ATTN:  Dr.  A.  Field 


46 


i 


t 


¥ 


T 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued* 


R 4-  D Associates 
ATTN:  Gerard  K 
ATTN: 

ATTN: 

ATTN: 

ATTN 
ATTN: 

ATTN. 

ATTN: 

ATTN: 


Schlcgel 
William  It.  Graham,  Jr. 
William  J.  Karzas 
Charles  Mo 
Leonard  Schlessinger 
S.  Clay  Rogers 
Robert  A . Poll 
Richard  It.  Schaefer 
Technical  Library 


The  Rand  Corporation 
ATTN:  Cullen  Cram 
ATTN:  Technical  Library 


Raytheon  Company 
ATTN:  Library 

ATTN:  Gajanan  II.  Joshi,  Radar  gvs. 


Lab. 


Raytheon  Company 

ATTN:  James  R.  Wcckback 
ATTN:  Technical  Library 


DEPARTMENT  OF  DEFENSE  ''ONTRACTOItS  (Continued) 


Inc. 


Science  Applications, 
Huntsville  Division 

ATTN:  Noel  It.  Byrn 
ATTN:  Technical  Library 


Science  Applications,  Inc. 
ATTN:  James  R.  Hill 
ATTN:  R.  Parkinson 
ATTN:  Richard  L.  Knight 


Simulation  Physics,  Inc. 
ATTN:  John  R.  Uglum 


Simulation  Physics,  Inc. 
ATTN:  It.  Little 


The  Singer  Company 

ATTN:  Irwin  Goldman,  Eng.  Management 
ATl’N:  Technical  Library 


Singer  Information  Systems  Network 

ATTN:  Technical  Information  Center 
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RCA  Corporation 

Government  4-  Commercial  Systems 
ATTN:  Technical  Library 

RCA  Corporation 

Government  4 Commercial  Systems 
ATTN:  Andrew  L.  Warren 
ATTN:  Technical  Library 

RCA  Corporation 
Camden  Complex 

ATTN:  E.  Van  Kcurcn,  13-5-2 
ATTN:  Technical  Library 

Rockwell  International  Corporation 
ATTN:  K.  F.  Hull 
ATTN:  Donald  J.  Stevens,  FA-70 
ATTN:  James  E.  Bell,  HA-10 
ATTN:  Technical  Library 

Rockwell  International  Corporation 
Space  Division 

ATTN:  TIC,  D/j  1-092,  AJ01 

Rockwell  International  Corporation 
ATTN:  T.  B.  Yates 

Sanders  Associates,  inc. 

ATTN:  Moe  L.  Aitcl,  NCA  1-3230 
ATTN:  Technical  Library 

Science  Applications,  Inc. 

ATTN:  Frederick  M.  Tesche 

Science  Applications,  Inc. 

ATTN:  William  L.  Chadscy 

Science  Applications,  Inc. 

ATTN:  Lewis  M.  Linson 
ATTN:  B.  H.  Fishbinc 
ATTN:  S.  J.  Dalich 
ATTN:  J.  N.  Wood 
ATTN:  R.  Fisher 
ATTN:  Technical  Library 


Southern  Research  Institute 
ATTN:  C.  Pears 

Sperry  Microwave  Electronics  Division 

Sperry  Rand  Corporation 

ATTN:  Technical  Library 

Sperry  Rand  Corporation 

Umvac  Division 

Defense  Systems  Division 

ATTN:  Technical  Library 

Sperry  Rand  Corporation 

Sperry  Division 

Soet.j  Gyroscope  Division 
ATTN:  Paul  Marraffino 
ATTN:  Charles  L.  Craig,  EV 
ATTN:  Technical  Library 

Sperry  Rand  Corporation 

Sperry  Flight  Systems  Division 
ATTN:  D.  J.  Keating 
ATTN:  Technical  Library 

Stanford  Research  Institute 
ATTN:  Philip  J.  Dolan 
ATTN:  William  C.  Taylor 
ATTN:  Arthur  Lee  Whitson 
ATTN:  George  Carpenter 
ATTN:  Dr.  G.  Abrahamson 
ATTN:  SRI,  Library,  Rm.  G-021 

Stanford  Research  Institute 

ATTN:  MacPhcrson  Morgan 
ATTN  Technical  Library 

Sundstrand  Corporation 

ATTN:  Curtis  B.  White 

Systems,  Science  4-  Software 
ATTN:  Technical  Library 

Systems,  Science  4-  Software,  Inc. 
ATTN:  Dr.  G.  Gurtman 
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DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

Systems,  Science  t Software,  Inc. 

ATTN:  Andrew  R.  Wilson 
ATTN:  Technical  Library 

Systron-Donner  Corporation 
ATTN:  Technical  Library 

Texas  Instruments,  Inc. 

ATTN:  Donald  J.  Manus,  MS  72 
ATTN:  Gary  F.  Hanson 
ATTN:  Technical  Library 

TRW  Semiconductors 

Division  of  TRW,  Inc. 

ATTN:  Technical  Library 

TRW  Systems  Group 

ATTN:  William  H.  Robinette,  Jr. 

ATTN:  Fred  N.  Holmquist,  MS  R1-202S 
ATTN:  Benjamin  Sussholtz 
ATTN:  Paul  Molmud,  R1-119G 
ATTN:  Aaron  H.  Narevsky,  IU-21 1 1 
ATTN:  A.  M.  Liebschutz,  R1-11G2 
ATTN:  Lillian  D.  Singletary,  Rl-1070 
ATTN:  Technical  Information  Center,  S-1930 
ATTN:  Robert  M.  Webb,  MS  Rl-1150 
ATTN:  Richard  H.  Kingsland,  Rl-2151 
ATTN:  Dr.  D.  Jortncr 

TRW  Systems  Group 

San  Bernardino  Operations 
ATTN:  John  E.  Dahnkc 
ATTN:  H.  S.  Jensen 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 

TRW  Systems  Group 

8 cy  ATTN-  Technical  Library 

United  Aircraft  Corporation 

Hamilton  Standard  Division 
ATTN:  Technical  Library- 

United  Technologies  Corporation 

Norden  Division 

ATTN:  Technical  Library- 

Victor  A.  J.  Van  Lint  (Consultant) 

ATTN:  V.  A.  J.  Van  Lint 

Varian  Associates 

ATTN;  Howard  R.  Jory,  A- 109 
ATTN:  Technical  Library 

Wcstinghousc  Electric  Corporation 

Astronuclear  Lalxiratory 

ATTN:  Technical  Library 

Wcstinghousc  Electric  Corporation 

Defense  & Electronic  Systems  Center 
ATTN:  Henry  P.  Kalapaca,  MS  3525 
ATTN:  Technical  Library 

Wcstinghousc  Electric  Corporation 

Research  1-  Development  Center 
ATTN:  Technical  Library- 
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